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^ " Abstract 

A novel interpretation of MOND is presented. For galactic data, in addition to Newtonian 
CN ■ acceleration, there is an attractive acceleration peaking at Milgrom's parameter clq. The 

J>-^ peak lies within experimental error where clq = cHq/2tt; Hq is the present-time value of 

the Hubble constant and c the velocity of light. The physical interpretation of this relation 
and its connection to Dark Energy are discussed. 






Oh 



X 



PACS numbers: 04.50. Kd, 98.62.Dm 



"p^; 1 Introduction 

It has long been known that rotation curves of galaxies disagree with Newton's law. Famaey 
and McGaugh have recently provided a review of all aspects of the data, with an exhaustive list 
of references [I] . The most relevant astrophysical data will be discussed here. 

In 1983, Milgrom proposed a modification of Newtonian mechanics called MOND (Modified 
£>v Newtonian Dynamics) [2], [3]. This has three distinct features. Firstly, the observed total 



acceleration a depends on the Newtonian acceleration g N as 

a = WMx); 



> 

H is an empirical smooth function of x — a / a o] a o is a universal constant ~ 1.2 x 10 10 ms 2 

?2 ■ for all galaxies. The second feature is that for vanishingly small g^, 

^j.' ■ a -> y/a g N . (2) 

o' 

Thirdly, galactic dynamics are invariant under scaling of time and space for very low g^\ 



(£,r)-> A(£,r), (3) 



where A is a scaling factor [3]. A star with rotational velocity v in equilibrium with centrifugal 
force satisfies 

2 j aoGM 

v l r =^—^' ( 4 ) 

G is the gravitational constant and M the galactic mass within radius r; so r cancels and 

v 4 = a GM. (5) 

This agrees well with the empirical Tully-Fisher relation between observed velocities at the 
edges of galaxies and their luminosities [5]. McGaugh showed that this relation applies over > 5 
decades of galactic masses from 10 6 to ~ 10 12 M Q after including the mass of gas and dust in 
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each galaxy J5J. It is important that in MOND there is only one free parameter oo, fitted to all 
galaxies once a particular form for fi(x) has been chosen. 

In the ACDM cosmological model, the parameter ao does not appear. There are three recent 
papers pointing out that a shift of paradigm is required to describe the data more precisely [TJ, 
[H]) IS]- Kroupa et al. and Famaey et al. discuss simulations based on the standard cosmological 
model ACDM. They both point out the prediction that the main Dark Matter halo hosting a 
large galaxy like the Milky Way should produce 100-600 roughly isotropic sub-halos. However, 
the Milky Way has only 24 satellites (and Andromeda ~ 28 [10]) which are highly correlated 
in both radial and momentum phase space. They lie mostly within rotationally supported thin 
discs. This has been strongly confirmed by Ibata et al. [11] . These authors all conclude that 
ACDM is at best incomplete and misses essential physics. There is a further interesting paper 
of Liighausen et al. concerning a polar ring galaxy [T2]. This galaxy has a small bright gas-poor 
disc with a large central bulge, but in addition an orthogonal gas-rich disc, referred to as a polar 
ring. They show that the observed velocities in both discs are well predicted by MOND. 

Another ACDM prediction is that dwarf galaxies formed from tidal material during galaxy 
encounters cannot contain substantial amounts of Dark Matter; MOND fits them well. McGaugh 
and Milgrom [S] make an important comment that tidal effects of large galaxies on their satellite 
galaxies must be taken into account in drawing conclusions about the satellites. The objective 
of the present paper is to search for an explanation of the observed discrepancies. 

The conventional ACDM paradigm is that Dark Matter condenses gravitationally and then 
galaxies form inside this condensate. This is a two-step process. The procedure adopted here is 
to use commonly used forms of Milgrom's fi function to determine the non- Newtonian component 
of the acceleration observed at the edges of galaxies; it peaks at or close to ao where it is bigger 
than gjv by a large factor ~ 171. This acceleration is then integrated to determine the associated 
energy function. The result fits naturally to a Fermi function with the same negative sign as 
that of gravity. The Fermi function lowers the total energy by 0.5 GM at radius r where g^ 
reaches ao; here M is the mass within radius r . This Fermi function takes the same form as 
the energy gap observed in doped semiconductors. It is adopted as the effective Hamiltonian 
for the non-Newtonian interaction. It is interpreted as evidence that a Fermi-Dirac condensate 
forms in the graviton-nucleon interaction (not the gravitational interaction itself) near radius 
r . There are then four further clues which fit like a glove to the existence of this condensate. 
One is that 2iia = cH within experimental errors; here c is the velocity of light and H is the 
local value of the Hubble acceleration. This results in a direct relation between ao and Dark 
Energy. Dark Matter is no longer needed in this one-step process. 

Most discussions of galactic rotation curves use the asymptotic form of the acceleration. Here 
attention is focussed on accelerations close to a , between ~ 10~ 8 and ~ 10~ 12 m s~ 2 . 

An independent observation of the role of ao has appeared in globular clusters within the last 
few years. These spherical clusters of stars have dimensions of a few parsec (pc), i.e. a factor 
~ 10 4 smaller than the Milky Way. They are believed to be the remnants of dwarf galaxies, 
some old and some quite young. Scarpa et al. reported initially on two globular clusters situated 
16-19 kpc from the Milky Way [13] . The equilibrium of such clusters is controlled by Jeans' 
Law. Scarpa et al. traced the velocity dispersion of 184 stars at large radius, identified as being 
members of one globular cluster (rather than interlopers), and 146 stars in the second cluster. 
The velocity dispersion is maximal at the centre of each cluster. They traced it to a radius twice 
that where gx reaches ao- Velocity dispersions deviate rather abruptly from g^ as it decreases 



through a . Tidal heating by the Milky Way varies as r -3 , and is at least one order of magnitude 
smaller, making its effect negligible. 

Scarpa et a!, have made observations of a further 6 globular clusters. Hernandez and Jimenez 
give the algebra relating velocity dispersions of stars to Newtonian acceleration using Jeans' Law 
[T5j. Hernandez, Jimenez and Allen report a detailed study of the velocity dispersion profiles 
of all 8 globular clusters [T5] . Like Scarpa et al., they conclude that tidal effects are significant 
only at radii larger by factors 2-10 than the radius where MOND flattens the curves. They 
also show that the velocity dispersion a varies with the mass M of the cluster as M -4 within 
errors; this is the expected analogue of the Tully-Fisher relation arising from Jeans' Law. This 
result is independent of luminosity measurements used in interpreting galactic rotation curves. 
In galaxies, the mass M within a particular radius is not easy to determine, and is usually taken 
as the mass where rotation curves flatten out. Further study of globular clusters is desirable. 

Differences between the \x functions used for MOND are illustrated in Fig. 19 of the review 
of Famaey and McGaugh pQ. The smoothest form, given by Milgrom [16], is used here: 

H(X) = V 1 + V(4X 2 ) - V(2x). (6) 

The algebraic form with which globular clusters are fitted is closely consistent with this equation. 

The layout of this paper is as follows. Section 2 presents in Fig. 1(a) the observed acceleration 
on the logarithmic scale of gw adopted by Milgrom [16] ; it gives the algebra for the chosen form 
of fJ>(x), where x = log 10 gAr. From this, the difference between the observed acceleration and 
(?tv is derived and shown in Fig. 1(b). This is then integrated to find its associated energy 
function W(x), Fig. 1(c). I follow the convention that the zero of the Newtonian potential is 
taken locally, ignoring the effects of Dark Energy over the Universe as a whole. The height of 
the Fermi function may be interpreted as an energy gap. 

There is further evidence of quantum mechanics at work. The asymptotic form of the total 
acceleration can likewise be integrated. In Section 3, this is shown to lead to a weak logarithmic 
tail to the Newtonian potential. The interpretation is that quantum mechanical mixing between 
the Newtonian potential and the asymptotic form allows the wavelength of gravitons trapped 
in the Newtonian potential to expand. This lowers the zero-point energy. 

These two indications of quantum effects should not be surprising. Hawking has demonstrated 
that Quantum Mechanics plays an important role in the physics of Black Holes [17], where 
the acceleration is very large. A standard Particle Physics result is that the graviton-nucleon 
interaction should obey a dispersion relation as a function of the amplitude of g^ (with a non- 
analytic logarithmic term in addition). Gravitons with a wavelength of galactic dimensions 
are at the extreme infra-red end of the spectrum. With such a wavelength, a single graviton 
interacts coherently with nucleons over complete clusters of stars. 

Section 3 examines the idea that there is quantum mechanical mixing between g^ and the 
Hubble acceleration. This mixing is a well known effect in Particle Spectroscopy for two eigen- 
functions having different basis states. The standard treatment is to use the Bogoliubov-Valatin 
transformation [T5] [TJJ], which leads to the Breit-Rabi equation for the mixing. 

Section 4 relates the observed asymptotic acceleration to Dark Energy and a subtle connection 
with the Hubble acceleration. Subsection 4.1 then draws an analogy with Condensed Matter 
Physics and the well understood Symmetry Breaking in Particle Physics. Section 5 makes 
extensive comments on avenues for further work and Section 6 summarises conclusion. 



1 THE MODEL NEAR a 

From equation (6), 



H = g N / a 

{g N /a + a /2a) 2 
g 2 N + a g N 




Oo 

4a 2 2a 
1 + (a /2a) 2 



(7) 



(9) 



The cancellation of the two terms a /2a is typical of a dispersion relation and emerges as an 
important point in Section 3. 
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-x = -log g 
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Figure 1: (a) The full curve shows the result of equation (10); the dashed line shows gw, and 
the dotted one a straight line given by ^a Q g N ; (b) the peak arising from curvature of the full 
curve in (a); (c) full curve: the energy derived from (b); the dashed curve is discussed in Section 
3; vertical dashed lines mark qq. 



Milgrom plots the log of total acceleration against the log of gravitational acceleration for 
systems of different masses [20]. It is inconvenient to use axes to the base 10, so they will be 
replaced in the algebra which follows by conversion to ln e . From Eq. (9) 



From equation (10), 





V = ln e yjgjf + a g N ; 




x = ln e g N . 


e 2 y 


= e 2x + a e x 


dy/dx 


= {e x + a /2){e x + a )- 1 


d 2 y/dx 2 


= (a /2)e x (e x + a y 2 


d 3 y/dx 3 


= (a /2)e x (a -e x )/(e x + a )- 3 ; 



(10) 

(11) 



(12) 
(13) 
(14) 
(15) 



d 3 y/dx 3 goes to zero at g — ao- [There are higher order terms too.] The curvature d 2 y/dx 2 has 
a maximum value at x = a , where d 2 y/dx 2 = 1/8. 



In galaxies, perturbations arise from thermal and pressure effects. Some narrow structures 
are observed in rotation curves, for example from bars in large galaxies. These may be fitted 
empirically using Poisson's equation for mass structures appearing in the gravitational potential. 
The MOND component 'rides' such structures smoothly, see Figs. 21 and 29 of Famaey and 
McGaugh pp. An immediate question is why Planck's constant does not appear in results for 
galactic rotation curves. The reason is that galaxies are noisy enough to hide it. 

Using ao = cHq/2tt, the maximum curvature is expected at ao = (1.113 ± 0.046) x 10 -10 
m s~ 2 , i.e. log 10 ao = —9.953. McGaugh summarises a large number of papers comparing the 
Tully-Fisher relation with models of galaxy formation [21]. He concludes that gas rich galaxies 
give the best determination of the baryonic masses of galaxies: a = (1.3 ± 0.3) x 10~ 10 m s~ 2 . 
A slightly lower value 1.22 ± 0.33 x 10~ 10 ms~ 2 is found by Gentile, Famaey and de Blok [22J. 
Fig. 1(b) shows d 2 y/dx 2 as a peak dW/dx in the acceleration at a , marked by the dashed line; 
here W is the energy of the 'extra' contribution. Results are insensitive to the precise value of 
ao, so figures and arithmetic are simplified by setting ao = 10~ 10 m s -2 . The full curve is well 
approximated by a Gaussian for the acceleration: 

d 2 y/dx 2 = 0.125 exp-[ 7 (x-a ) 2 ] (16) 

with 7 = 1.175, i.e. the Gaussian drops to half-height at 9.6% of the value of x at the peak in 
Fig. 1(b). The conclusion is that galaxies have considerable stability. 

2.1 Alternative forms for \i 

Other forms for fi(x) have been used to fit galactic rotation curves. Two common examples 
are (A) /z(x) — x/(\ + X)i (B) x/\A + X 1 ■ Both forms have the effect of moving the peak of 
the curvature slightly: from x = —10 to —9.7 for A and to —9.85 for B. The height of the 
peak for A is scaled by a factor 0.78 and the curve becomes correspondingly wider, resulting in 
a tail reaching 0.012 at x = —8.2 and —11.8. The result for Fig. 1(c) is that the top of the 
Fermi function is 4% of GM lower, and the bottom of the curve higher by the same amount, 
but the central part of the Fermi function is unchanged and still centred very close to x — —10. 
Form B produces the converse effect: a higher, narrower peak in Fig. 1(b) and a Fermi function 
beginning closer to the top and finishing closer to the bottom of Fig. 1(c), but with its central 
section undisturbed. 

A detail is that one might consider the possibility that the smooth curve of Fig. 1(a) could 
be derived from the dashed and dotted lines. This would give rise to a sharp cusp in Fig. 1(b). 
Such cusps are known in Particle Physics, but arise only at the opening of phase space for a 
new reaction channel [23]. However, no such threshold exists in galactic phenomena. 

In Fig. 1(a), the dashed line corresponding to ajv intersects at x = — 10 with the straight 
line fitted to y/aogN- This suggests that there is a cross-over between two eigenf unctions corre- 
sponding to two asymptotic regimes. This is a familiar result in Particle Physics. Historically, it 
was first observed in variations of atomic energy levels in a magnetic field. It arises there from 
differences in magnetic moments for different levels, some of which share one component. The 
result for that case is sketched in Fig. 2. The two dashed lines show the variation of energy 
eigenstates with magnetic field. 




Figure 2: Sketch of two crossing atomic lines; full lines show them including mixing, dashed 
lines without; E\ and E 2 label the convention for eigenvalues in the case of no mixing. 



For the two mixed states, 



#22*2 + ^1 = E^ 2 - 

where V is the mixing Hamiltonian. The solution of these coupled equations is 



(H 11 -E)(H 22 -E)-V 2 



0. 



(17) 
(18) 



(19) 



This equation was first derived in 1931 by Breit and Rabi [21]. The same formalism describes 
mixing between the three neutrinos u e , v^ and v T and also the CKM matrix of QCD. For galaxies, 
classical expectation values < H\\ > and < H 2 2 > are to be substituted into Eq. (19). 




Figure 3: Axes x, y, x' and y' '. 



The relevance to galactic physics is that Fig. 1(a) includes mixing between Newtonian gravity 
(to the left of the figure) and the Hubble acceleration (to the right). In order to reproduce the 
symmetric form of Fig. 2, the Bogoliubov transformation is needed. It rotates Fig. 1(a) by 
35.78° anti-clockwise; this is the mean angle of the dashed and dotted lines with respect to the 
x-axis. The rotation is illustrated in Fig. 3; it is about the point x = —10, y = —10, where the 
two straight lines of Fig. 1(a) cross: 



x' = (a; + 10) cos /3- (y + 10) sin/3 
y ' = ( x + 10) sin ft + (y + 10) cos (3. 



(20) 
(21) 



Substituting Eq. (10) gives an exact expression for the curve in x',y' axes. 

It is also convenient to re-express gjy and the Newtonian energy directly in terms of x; from 
Fig. 1(a) 

H u = E l = -GM/r = -V~GMe x/2 (22) 

H 22 = E 2 = VGMe(x) (23) 

V = VGMW(x'), (24) 

where e refers to an energy possibly given by the Hubble acceleration; however, in practice, the 
effect of the Hubble acceleration over the radius of the Milky Way is < 2 x 1CT 4 , which can 
be neglected. Note that the factor \/GM appears in all three equations, so as to conform with 
Milgrom's scaling law. The common approximation is also used that the gravitational energy 
of a disc galaxy is —GM/r, where M is the mass inside radius r. This is accurate to < 1% at 
the large values of r of interest. The two solutions of the Breit-Rabi equation are 



E = ^±^ ± \(^r^-) 2 + V\ (25) 



3 The relation to a Fermi function 

The acceleration differs in x and x' axes, but the scalar quantity W is independent of axes. The 
standard form of the Fermi function is 



W(x) oc 



l+e ^( E ji L , 



(26) 



where Ep is the energy at the centre of the Fermi function and (3 is a fitted constant. This 
refers to a situation where there is a discrete energy gap, e.g. in a perfect semi-conductor or 
superconductor. W(x) can be obtained by numerical integration of y from Eq. (10) and is 
shown as the full curve in Fig. 1(c). 

3.1 A long-range logarithmic tail to the Newtonian potential 

If there were no quantum mechanical mixing between Newtonian acceleration and the conden- 
sation mechanism, there would be no structure at the crossing point and, more serious, no 
explanation for the term y/a gN- 

Asymptotically, the total acceleration, taken from MOND, is 



a = a ^g N /a . (27) 

Since g N = GM/r 2 , 



a -+ ^/GMa /r. (21 

Taking this as —d(j)/dr, where is a potential induced by the mixing, 



-^-JGMaoHr/n). (29) 



Here r\ is the mean radius for this term. The value of r\ is necessarily very close to the much 
larger dip caused by dW/dx. Because a® ~ 10 -10 , is very small. However, it does explain the 
asymptotic straight-line at the right-hand edge of Fig. 1(a) and contributes to e{x) in H 2 2- 

The interpretation of this term is that mixing between the Newtonian potential and the con- 
densation mechanism expands the wave-length of gravitons trapped in the Newtonian potential 
and lowers their zero-point energy. An analogy is the covalent bond in chemistry. In the hy- 
drogen molecule, each of two electrons is attracted to two protons (which themselves repel one 
another). This increases the wave-length of the electrons and reduces their zero-point energy. 

Sobouti noticed this long-range logarithmic component and wrote two closely related papers 
[25] . |26j . These included effects of General Relativity. This is interesting but not strictly 
necessary at present. The problem was solved by Sobouti as a power series with additional 
empirical terms proportional to 1/r and 1/r 2 . These terms are now replaced by our equations. 
There is experimental evidence that the Milky Way has a halo extending to ~ 100 kpc [27]. 
However, this could be due to gas and dust shared with the local cluster of galaxies. So this 
halo is presently ambiguous. 

3.2 Solving the Breit-Rabi equation 

The binding energy W is defined with the same sign as E-\_. From equations (25) and (24), 



IE = E 1 + e(x)±^(E 1 -e(x)) 2 + AV 2 (30) 

Vgm 



ln P 10 



(-e x/2 + e(x) ± ^J(e x / 2 -e(x)) 2 + 4wA (31) 

2 dE = VGM ( Q 5e , /2 | ± (e*/ 2 - e(x))(0.5e-/ 2 - de(x)/dx) + 4WdW/dx \ 

dx ln e 10 ^ ' ' ' ~'~ yf(e*/ 2 -e(x)) 2 + 4W 2 / 

The factor ln e 10 allows for the fact that Fig. 1(a) has been drawn using axes which use logarithms 
to the base 10. For the upper branch of the solution, the minus sign for the term involving the 
square root is required to reproduce the usual Newtonian potential. The equations simplify 
considerably if the very small term e(x) is neglected. 

A Bose-Einstein condensate does not fit the data. In this case, dW/dx should vary as T 3 ' 2 at 
the peak of the acceleration [28]. For positive x near x = 0, the relation of the energy function 
to kT gives, 

W(x) = -5(l + |x| 3/2 exp-7'x 2 ) (33) 

dW(x)/dx' = B(1.5\x\ - 5 -27'|x| 5/2 )exp-7'x 2 ; (34) 

for negative x, the opposite sign of |x| is needed in W(x). The second term in dW{x)/dx wrecks 
the x dependence, which fails to fit the observed peak. Near x = it has a square root variation 
with x and then, when the second term of Eq. (34) overtakes it, the curve turns downwards. 
This rules out a Bose-Einstein condensate. 

The depth of the Fermi function is -0.5GM. The magnitude of this term can be traced to 
the factor 2 difference in slope of gjy and that of the asymptotic form ^a g N . 

Let us now return to Figs. 1(b) and (c). Here there is a slight complication. Fig. 1(b) is 
the acceleration measured in x,y axes. However, the rotation to x',y' axes requires that 7 of 
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Eq. (16) is increased to 1.852. In addition, there is a small visible displacement of the centre 
of curvature in Fig. 1(a) by an offset of —0.203 x 10~ 10 in x. What then emerges from the 
Breit-Rabi equation is that g^ is rather small near x = a^ compared with that originating from 
the extra acceleration dW/dx'. This second term dominates by a large factor ~ 171 at the centre 
of the curve. This ratio falls by 50% at x = —10.6, to 30 at x = —11, then -6.0 at x — —11.5 
and ~ 1 at x = —12. Results for the 'extra' acceleration are symmetric about a except for the 
term ^a g^ of Eq. (10). The conclusion is that the curved part of Fig. 1(a) is the dominant 
feature near x = a where Newtonian gravitation is a rather small perturbation. 

Consider the effect of this result near the centre of the Fermi function at E — Ep in Fig. 
1(c). If we retain only the dominant terms in W and dW/dx, Eqs. (31) and (32) give 



dE/dx -> i- IdWldx (35) 

m e 10 

E -+ ^2W. (36) 



Apart from the factor y/GM/ln e 10, which is a normalisation factor, dE/dx may be interpreted 
as the modulus of a Breit-Wigner resonance with x-dependent width: 

BW = r W /2 

E-E F - iT(x)/2 K ' 

The energy W starts at zero because of local gauge invariance, and its central value is shifted 
downward by 0.25 GM. This accounts for the form of Eq. (8) of Section 2, where two terms 
(a /2a) 2 cancel. 

How can this effect be understood? A possibility is that the graviton acquires a small effective 
mass near oo- The fitted change to the gravitational acceleration is close to a Gaussian, as in 
Fig. 1(b). In subsection 2.1, the dependence of the fit on alternative forms of Milgrom's fi 
function was tested. Although acceleration curves change significantly, the Fermi function is 
affected only at the ends of the range x = 8 to 12 by at most ±4%. 

The conclusion from these results is that the central part of the Fermi function is stable, 
but can be perturbed at the edges. In superconductors, a coherence length was introduced by 
Pippard to account for the effects of defects beyond an experimentally observed range [29]. It 
appears that galaxies behave similarly. If gravitons acquire an effective mass, it appears at first 
sight that this will weaken gravity. However, remember that the wave-lengths of gravitons on a 
galactic scale are very large. The spectrum of gravitons near the edge of the galaxy is a convolu- 
tion of gravitons from the rest of the galaxy. Gravitons from the galactic centre become almost 
plane waves which can interact coherently over a large volume. They function like phonons 
coupling to Cooper pairs in a superconductor producing an energy gap. A coherence length like 
that introduced by Pippard can arise in many ways. Supernovae act as major perturbations, 
heating sizable volumes. It is also known that so-called chimneys and wormholes provide chan- 
nels through which currents of dust and gas flow. Furthermore, one of the remarkable features 
of galaxies is that they only grow to a certain size. The largest have masses of order 10 12 M Q . 
This can be attributed to the rapid fall-off of the 'extra' acceleration due to MOND close to 
x = —8 and —12. All of these effects point to a coherence length due to the variation of the 
Fermi function in Fig. 1(c). 



Let us return to a simpler issue, the missing lower branch of the Breit-Rabi equation. On this 
branch both W and dW/dx' change sign. The change of sign requires that this branch describes 
an excited state rather than a condensate. (Remember that energies of both gravity and W(x) 
are negative). Such an excited state is likely to decay on a time scale much less than that of 
galaxies, so it is unlikely that this branch will be observable. For those wishing to investigate 
this branch, the procedure is (a) to fit the upper branch as a function of x, (b) rotate to x',y' 
axes, (c) reverse the sign of y' to reach the lower branch, and (d) rotate axes back again to x, y. 
The best place to search for this branch is near the crossing point of Fig. 1(a). 

4 The relation to Dark Energy 

Experiment tells us that in galaxies, the asymptotic form of the acceleration is a/oo^/v- This 
leads to the question: what governs the asymptotic acceleration? 

If MOND successfully models the formation of galaxies and globular clusters, it raises the 
question of how to interpret Dark Energy. In a de Sitter universe, the Friedmann-Robertson- 
Walker model smoothes out structures using a ACDM function which models the gross features. 
These change over the lifetime of the Universe. However, if quantum mechanics governs indi- 
vidual galaxies, there will instead be fine structure in Dark Energy. It is logical that steps like 
Figs. 1(c) do not just average out, but instead accumulate over all of this fine structure. 

The standard treatment of Dark Energy is reviewed by Sami [30]. The acceleration falls 
initially due to Newtonian cosmology. However, at recent times, the acceleration increases. 
This acceleration can be explained naturally by the sum total of the fine structure over all the 
galaxies. The total acceleration is parametrised via the assumed time dependence of the metric 
on the Hubble acceleration. A similar suggestion along these lines has been advanced by Zhang 
and Li [31] using ideas based on entropic arguments. The 'present-day' Hubble acceleration is 
the local value and varies over the Universe according to the parametrisation by Dark Energy. 

There is a further argument pointing towards the idea that local fine structure is cumulative. 
Peebles and Nusser argue that galaxies condense more rapidly than the standard ACDM model 
predicts [32] . In particular, they point out that the Local Void contains far fewer galaxies than 
ACDM predicts statistically, while there is an unexpected presence of large galaxies on the 
outskirts of the Local Void. Their Fig. 1 is very persuasive in this respect. Only 3 galaxies 
are observed in the Local Void compared with 19 predicted. The Poisson probability for this 
result is of the order 10~ 5 from ACDM. Peebles and Nusser conclude: 'In short, the general 
sensitivity of galaxies to their environment is not expected in standard ideas. It would help if 
galaxies were more rapidly assembled so that they could then evolve as more nearly isolated 
island universes.' Later Peebles considered an additional empirical term added to the ACDM 
model, but commented that the change requires that Cold Dark Matter is cored rather than the 
expected cusped behaviour [33] . 

A natural explanation is that the Local Void gives no contribution to the Hubble mechanism, 
except for its three galaxies. The total energy E is then higher there. On the periphery of 
the Void, there is gas and dust which can form galaxies. This gas runs down the energy gap 
to enlarge galaxies forming there. Galaxies then grow by accumulation of dust and gas. The 
condensate at the edge of the galaxy acts as a funnel to collect gas and dust. Towards the 
centres of galaxies, they can develop in different ways depending on the angular momentum 
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L of the galaxy. For large L, they naturally form flat discs. For lower L, bulges develop and 
spiral arms; these features can vary according to whether the galaxy is supported by rotation or 
pressure or both. For really low L, large elliptical galaxies and dwarf spherical galaxies develop. 
If L is close to 0, galaxies collapse as quasars. It appears that AC DM is presently fitting all 
these different morphologies with a very flexible ansatz for Dark Matter. 

So far, Dark Matter has escaped experimental detection. There are speculations that it 
may take the form of sterile neutrinos with masses in the electron-volt range. At a recent 
conference on neutrino physics, Altarelli commented on this question j3l] . He argued that more 
than one sterile neutrino is disfavoured by stringent bounds arising from nucleo-synthesis. Also 
there is tension between LSND, MiniBoone and KARMEN experiments leaving little room for a 
signal. Rubbia et al. propose a new neutrino detector ICARUS-NESSIE to run at CERN with 
parameters optimised for finding sterile neutrinos [55] . 

There is also speculation about Grand Unified Theories in which massive right-handed neu- 
trinos couple weakly to left-handed light neutrinos. Such heavy neutrinos would have formed 
in the Big Bang before the generation of the Cosmic Microwave Background. If these heavy 
neutrinos survive and mix with light neutrinos, their lifetimes must be larger than the age of 
the Universe, otherwise there would be effects visible in Dark Energy. 

4.1 Other Condensates 

Well known examples of condensates are ferromagnets and anti-ferromagnets. In ferromagnets, 
spins align parallel with one another below the Curie point; in the absence of a magnetic field, 
the overall spin can lie in any direction. In anti-ferromagnets they spontaneously align anti- 
parallel. The ground-states of these systems have a lower symmetry than the Hamiltonian. 
This is a purely quantum effect, but difficult to calculate from first principles. However, neutron 
scattering experiments establish the structures. In astrophysics, the situation is difficult because 
one can only watch how galaxies evolve. 

In Particle Physics, we now know that the Strong Interaction is mediated by massless gluons 
obeying Chiral Symmetry, i.e. they do not discriminate between left-hands and right-hands. 
Below a mass of ~ 1 GeV, the gluon acquires an effective mass from its interaction with light 
quarks, which themselves have masses of ~ 4 and 9 MeV. The Electroweak Theory is constructed 
from a mixing between Electromagnetism and Weak Interactions carried by W and Z particles 
and probably the Higgs boson. A side-effect of this idea is that Chiral Symmetry is broken in 
the strong interactions for spin particles below 1 GeV. This idea was introduced by Gasser 
and Leutwyler [36] and today gives many quantitatively accurate results in meson spectroscopy. 
An important step in understanding the precise mechanism was made by Bicudo and Ribiero 
|37j . Their work finds the Bogoliubov-Valatin transformation to be an essential element. Above 
1 GeV, there is a cross-over in which Chiral Symmetry is largely restored and the quark model 
reigns supreme, though with small amounts of mixing with meson-meson and/or qqqq basis 
states. A recent paper of Pennington and Wilson gives details including figures showing the 
cross-over at 1 GeV [38] . 

A precise set of equations describing Chiral Symmetry Breaking is given by Cherney et al. [3H] ■ 
This work shows that gauge invariance requires that quarks must be treated as 'dressed' fermions, 
rather than bare fermions. Their conclusion is that a Yukawa potential appears explicitly in the 
qq interaction, their Eq. (61). This increases the fermion mass from a few MeV to the so-called 
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constituent mass, ~ 320 MeV. In their Eq. (72), a logarithmic term appears from Feynman 
diagrams where one pair of particles rescatter to themselves through a closed loop. The origin 
of Chiral Symmetry Breaking is that quarks have small masses; if these were zero, the pion 
would be massless. 

In summary, it is clear that there are three gauge fields: (i) gluons, (ii) electromagnetism 
linked to weak interactions, and (iii) gravity. Effective masses do arise in the first two. There 
has been a paper by Van Dam and Veltman [ID] claiming a no-go theorem preventing the 
graviton developing a mass. This theorem appears to refer purely to an isolated graviton obeying 
General Relativity. That is different to the present case where the condensate is a property of 
the graviton-nucleon interaction. 

There is an important point if gravitons develop an effective mass. In Particle Physics, there 
are relations between scattering processes such as ii + p —¥ ir + p, ?c~p — >■ tc~p and tt + tc~ — y pp. 
Amplitudes depend on E 2 — p 2 , where E and p are momenta in each reaction. These are 
denoted s, u and ^-channels by a convention introduced by Mandelstam. For massless gravitons, 
kinematics of the s-, t- and u- channels all meet at a point. However, if the graviton acquires an 
effective mass, it moves these three channels apart and changes the graviton-nucleon interaction. 
In particular, it affects the nucleon pole term in the w-channel, hence the structure of the nucleon. 
What is needed, but beyond reach at present, is to solve the Schwinger-Dyson equations for this 
case following the procedure of Pennington and Wilson for the im system. The obstacle is that 
the structure of the nucleon is not yet known precisely enough for a realistic calculation from first 
principles. Thomas developed a model in 1983 explaining features of deep inelastic scattering 
(i.e. large momentum and energy transfers in electron scattering from nucleons) which provided 
a qualitative explanation of the data in terms of a chiral quark model [H] . The idea is that a 
nucleon has a component in its wave function with a pion circulating round it with one unit of 
orbital angular momentum: iV — y [N + tt]l=i- This has been followed up in a recent paper of 
Burkhardt et al. which improves this work in specific ways [12] , but the final set of parameters 
is not yet available. 

Cosmology needs to be treated equally empirically. My suggestion is that Dark Energy is 
symptomatic of spontaneous symmetry breaking of gravitation to a de Sitter universe governed 
by the space and time components of the metric. The de Sitter universe has a lower symmetry 
than General Relativity. In galaxies, ao is the order parameter of a Fermi-Dirac condensate and 
is directly related to cH /2tt. The origin of the factor 2-7T is the same as that between Planck's 
constant and h. In both cases, what is involved is a transformation from spherical or cylindrical 
coordinates to Cartesian coordinates. 

5 Further work which is needed 

The model proposed here is precise and open to experimental test. The most important and 
simplest is that if Dark Matter is replaced by this model, it is obviously necessary to redo the 
parametrisation of Dark Energy so that it reproduces smoothly what has been parametrised 
as Dark Matter up to now. This does not necessarily require major modifications. The main 
point is to fit the third peak in the Cosmic Microwave Background. This requires cooperation of 
groups with the latest Planck data and techniques at their finger-tips. Despite the clues pointed 
out here that a Fermi-Dirac condensate explains galactic rotation curves, this could fail. It 
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would not be surprising if minor modifications are required in the refit to Dark Energy and new 
clues might emerge. It is necessary to include the logarithmic tail of the Newtonian potential 
found here. It is also essential to treat Voids so that they only contribute to Dark Energy via 
the galaxies observed there. Those may not be well resolved if they are very distant from us. It 
is also essential to fit the Hubble parameter and magnitudes of Typela supernovae as a function 
of red-shift z. 

Whatever develops at the edges of galaxies will affect both Strong and Weak Gravitational 
Lensing. However, unless data are used where a distant quasar transmits light through the 
periphery of a galaxy, where the effect of W(x) is large, the lensing effect will be rather small. 
The weak logarithmic tail of Newtonian gravitation also needs to be taken into account. 

What happens in clusters of galaxies needs detailed, laborious calculations using the formulae 
given here; since there are presently claims that Milgrom's formula does not correctly reproduce 
what happens in clusters, this could be revealing. A multi-body interaction between the Hubble 
mechanism and several galaxies is required. Quite apart from the attraction to acceleration ao, 
there are also strong tidal effects of the variety discussed by Kroupa |43] . 

The same remarks apply to the Bullet Cluster, which needs to be refitted with the model 
proposed here. The calculation becomes a two-centre problem. The Hubble acceleration couples 
to both galaxies, but they also couple to one another, modifying the zero-point energy; individual 
stars in galaxies communicate with one another not only through the Newtonian potential but 
via their Fermi functions. In the Bullet Cluster, each galaxy behaves in this way, but there 
will be coupling between stars "belonging" to each individual galaxy. There may be complex 
interactions between the two galaxies including resonance effects. There are new data on the 
Bullet Cluster jH] showing that there are many dwarf satellite galaxies in the cluster. 

6 Concluding Remarks 

The equations given here allow very little freedom - just the 4% perturbations allowed in the 
Fermi function at the top and bottom of Fig. 1(c). There are five clues which point to galaxies 
and globular clusters being quantum mechanical condensates. 

1. The phenomena of Fig. 1 appear on a logarithmic scale of g^. This is the form expected 
for the Partition function of Statistical Mechanics: 

Z = n *rr^ (38) 

where E s are energies of levels in a box and n denotes the numbers of quantum levels in 
the box. For a refresher course, see Schrodinger's clear exposition 



2. Fig. 1(a) can be interpreted in terms of quantum mechanical mixing between two crossing 
eigenstates; the rotation of axes accomodated by Fig. 3 is just what is expected from the 
Bogoliubov-Valatin transformation, a quantum effect. 

3. Fig. 1(c) is fitted naturally by a Fermi function with an energy gap 0.5 GM. 

4. The asymptotic form of the acceleration in Fig. 1(a) generates a logarithmic tail; this 
requires a quantum mechanical explanation. A Fermi-Dirac condensate fits the data; a 
Bose-Einstein condensate does not. 

13 



5. Long wave-length gravitons can explain the amplification of the amplitude forming the 
condensate. 

The lower branch of the Breit-Rabi equation has the opposite curvature to the upper one. 
The natural interpretation of this result is that it corresponds to an excited state which will 
decay rapidly and will not therefore appear in galactic rotation curves. 

The interpretation of the curve of Fig. 1(b) as an energy-dependent Breit-Wigner pole is that 
the graviton acquires an effective mass in the vicinity of a . How this originates is speculative 
but does not reflect on the five clues listed above; it depends on the nucleon structure function 
which is not yet well defined. 

The asymptotic form of the acceleration is y/a Q g. This is clearly associated with the Hubble 
acceleration and is related to Dark Energy. The standard approach to the Universe is the 
Friedmann- Robert son- Walker model where components of the metric are smoothed over local 
structures and appear in the radial and time components of the metric of General Relativity. 
My suggestion is that galaxies create fine-structure and the Friedmann- Robertson- Walker model 
should include into Dark Energy the sum over these structures. This may explain the agreement 
between a Q and cH /2n in our part of the Universe as well as its late-time acceleration. 

Those are the essential points. They do reflect successfully the facts which M OND parametrises. 
However, only the peripheries of galaxies are considered. The morphology at the centres of 
galaxies is a separate issue. 

There are two further small comments. Firstly, there have been many papers on the derivation 
of Newton's law from Entropic arguments. Many assume that one can assign a temperature T 
to the Hubble acceleration and a temperature 7\ for the gravitational potential by adding them 
as the sum of squares, i.e. the random phase approximation. This gives the result for the total 
acceleration: 



\9N + 9H) 2 -g 2 H (39) 

where gn is the Hubble acceleration H^r, and r is the distance to the centre of the gravitating 
object. This is a quite different approach to the one proposed here and leads to a very small 
effect. For the Milky Way, it is a maximum of 2 x 1Q~ a cHq/2tt. 

Secondly, Iorio has calculated that effects of MOND on the perihelia of planets in the Solar 
system are about a factor 10 below present experimental errors [16]. The logarithmic term 
arising from MOND will be of order a Q , which is very small and will vary exceedingly slowly 
over the solar system. It has been pointed out by Galianni et al. that it may be feasible to 
detect the MOND effect at acceleration aQ in the solar system near Lagrangian points where 
the accelerations of the Sun, Earth and Moon cancel out |47) . If measurements of sufficient 
accuracy could be made, they might confirm or modify MOND as a model of the behaviour 
of galaxies; secondly, they have the important potential to measure the shape of the response 
function through the region where the bend appears in Fig. 1. Further study of globular clusters 
would also provide information on the rotation curves of galaxies. 
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